function g = constrains(x)
    g = zeros(4,1);
    g(1) = - 11 + x(1)^2 - 6 * x(1) + 4 * x(2);%换成小于等于
    g(2) = - x(1) * x(2) + 3 * x(2) + exp(x(1) - 3) - 1;
    g(3) = - x(1);
    g(4) = - x(2);
end
%约束条件